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We study in cranked Nilsson plus random phase approximation shape transitions in fast rotating 
nuclei undergoing backbending, more specifically ^®®Dy and ^®^Yb. We found that a backbending in 
^®®Dy is correlated with the disappearance of the collective, positive signature 7 -vibrational mode in 
the rotating frame, and, a shape transition (axial—> nonaxial) is accompanied with a large acquiring 
of the 7 deformation. We show that such a shape transition can be considered as a phase transition 
of the first order. In '^^^Yb the quasiparticle alignment dominates in the backbending and the shape 
transition (axial^ nonaxial) is accompanied with a smooth transition from zero to nonzero values 
of the 7 deformation. We extend the classical Landau theory for rotating nuclei and show that the 
backbending in ^^^Yb is identified with the second order phase transition. A description of spectral 
and decay properties of the yrast states and low-lying excitations demonstrates a good agreement 
between our results and experimental data. 

PACS numbers: 21.10.Re,21.60. Jz,27.70.-|-q 


I. INTRODUCTION 

Backbending is a paradigm of structural changes in a nucleus under rotation. A sudden increase of a nuclear 
moment of inertia in yrast rotational band at some critical angular momentum or rotational frequency, discovered few 
decades ago Q, continues to attract a considerable attention. There is a general persuasion that this phenomenon 
is a result of the rotational alignment of angular momenta of a nucleon pair occupying a high-j intruder orbital near 
the Fermi surface (see textbooks ill and references therein). However, it is understood just recently, that this 
point of view may obscure different mechanisms if applied to nuclei with relatively small axial deformation at zero 
spin. In particular, Regan et al 0 proposed that the backbending observed in a number Cd, Pd and Ru nuclei can 
be interpreted as a change from vibrational to rotational structure. In our preliminary report we proposed to 
consider the backbending in TV ~ 90 as a result of the transition from axially symmetric to nonaxial shape due to a 
disappearance of a collective 7 -vibrational mode in the rotating frame. 

The analogy of the backbending phenomenon with a behaviour of superconductors in magnetic field, noticed in 
Refi, is rising a desire to apply Landau theory of phase transitions Q to nuclei. We recall that Landau’s theory deals 
with second order phase transitions, when different macroscopic phases become indistinguishable at the transition 
point. Whether it takes place for rotating nuclei is an open question. Nuclei are finite systems and phase transitions 
should be washed out by quantum fluctuations. Nevertheless, long ago Thouless 0 proposed to distinguish two kinds 
of ’’phase transitions” even for nuclei. Such phase transitions may be connected with shape transitions, for example, 
from spherical to deformed or axially deformed to nonaxially deformed shapes. This idea was put ahead in the analysis 
of shape transitions in hot rotating nuclei In this case the statistical treatment of the finite-temperature mean 
field description had provided a justification for an application of the Landau theory for nuclei. Within this approach, 
a simple rules for different shape-phase transitions were found as a function of angular momentum and temperature. 

Recently, quantum phase transitions, that occur at zero temperature as a function of some nonthermal control 
parameter, attract a considerable attention in various branches of many-body physics, starting from low-dimensional 
systems ^ to atomic nuclei and molecules Until the present a major activity in the study of shape-phase 

transitions for nuclei in the ground state at zero temperature is carried out within the interacting boson model (IBM) 
(see, for example, R,efsll2L lllll and references therein). The model naturally incorporates different symmetry limits 
associated with specific nuclear properties 0 . While the IBM can be easily extended to a thermodynamical limit 
TV —> 00 , which is well suitable for the study of phase transitions, the analysis is rather oversimplified. For example, 
the model does not take into account the interplay between single-particle and collective degrees of freedom in even- 
even nuclei. A general trend found for the ground shape transitions is less affected by this interplay. However, it may 
be crucial for the study of quantum phase transitions in rotating nuclei, where statical and dynamical properties are 
coupled. As we will see below, the above interplay determines the type of quantum shape-phase transitions in rotating 
nuclei. It elucidates also the behaviour of low-lying excitations, specifically related to the shape transitions at high 
spins. Among such modes are 7 -excitations and wobbling excitations, which are related to the nonaxial shapes. 
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The nuclear shell model (SM) treats the single-particle (s.p.) and collective degrees of freedom equally and appears 
to be extremely successful in the calculation of the backbending curve in light nuclei However, the drastic increase 
of the configuration space for medium and heavy systems makes the shell model calculations impossible. In addition, 
one needs some model consideration to interpret the SM results. On the other hand, various cranking Hartree-Fock- 
Bogoliubov (HFB) calculations (cf [la US Ell) provide a reliable analysis of the backbending for medium and heavy 
systems. As a rule, low-lying rotational bands are described within the cranking model with a principal axis (PAG) 
rotation. For the PAG rotation each single-particle (quasiparticle) configuration corresponds to a band of a given 
parity and a signature E3- In the HFB calculations the backbending is explained as a crossing of two quasiparticle 
conhgurations with different mean field characteristics. 

It is well known, however, that a mean field description of finite Fermi systems could break spontaneously one 
of the symmetries of the exact Hamiltonian, the so-called spontaneously symmetry breaking (SSB) phenomenon 
(see RefsIH 113 for a recent review on the SSB effects in rotating nuclei). Obviously, for finite systems quantum 
fluctuations, beyond the mean held approach, are quite important. The random phase approximation (RPA) being an 
efficient tool to study these quantum fluctuations (vibrational and rotational excitations) provides also a consistent 
way to treat broken symmetries. Moreover, it separates collective excitations associated with each broken symmetry 
as a spurious RPA mode and hxes the corresponding inertial parameter. This was recently demonstrated in fully 
self-consistent, unrestricted Hartree-Fock (HF) calculations for other mesoscopic system, two-dimensional quantum 
dots with small number of electrons and a parabolic conhnement 1^ . For large enough values of the Goulomb 
interaction-conhnement ratio Rw the HF mean held breaks circular symmetry; the electrons being localized in specihc 
geometric distributions. Applying next the RPA it was shown that the broken symmetry corresponds to appearance 
of spurious (with a zero energy) RPA mode (Nambu-Goldstone mode). And this mode can be associated with a 
rotational collective motion of this specihc (deformed) electron conhguration (see and references therein), 

which is separated from the vibrational excitations. Thus, a self-consistent mean held calculations combined with 
the RPA analysis could be useful to reveal structural changes in a mesoscopic system, i.e., to detect a quantum 
shape-phase transition. 

In contrast to quantum dots, in nuclei, a nucleon-nucleon interaction is less known. Mean held calculations with 
effective density dependent nuclear interactions such as Gogny or Skyrme forces or a relativistic mean held approach 
still do not provide sufficiently accurate single-particle spectra to obtain a reliable description of experimental char¬ 
acteristics of low-lying states (cf jl^l). The RPA analysis based on such mean held solutions is focused only on the 
description of various giant resonances in nonrotating nuclei, when the accuracy of single-particle spectra near the 
Fermi level is not important (cf (l^l). Furthermore, a practical application of the RPA for the nonseparable effective 
forces in rotating nuclei requires too large conhguration space and is not available yet. A self-consistent mean held 
obtained with the aid of phenomenological cranked Nilsson or Saxon-Woods potentials and pairing forces is quite 
competitive up to now, from the above point of view. These potentials allow to construct also a self-consistent resid¬ 
ual interaction neglected at the mean held level. The RPA with a separable multipole-multipole interaction based 
on these phenomenological potentials is an ehective tool to study low-lying collective excitations at high spins (cf 

mill). 

It was demonstrated recently, in an exactly solvable cranking harmonic oscillator model with a self-consistent 
separable residual interaction 1^ , that a direct correspondence between the SSB effects of the rotating mean 
held and zero RPA modes can be established in a rotating frame if and only if mean held minima are found self- 
consistently. Thus, it is self-evident that the analysis of the SSB ehects and RPA excitations for realistic potentials 
requires a maximal accuracy of the fulhllment of self-consistency conditions. In Refi, we proposed a practical method 
to solve almost self-consistently the mean held problem for the cranked Nilsson model with the pairing forces in order 
to study quadrupole excitations in the RPA. In the present paper we discuss all the details of our method and analyse 
the backbending in ^®®Dy and ^®^Yb. We thoroughly investigate the positive signature quadrupole excited bands as 
a function of the angular rotational frequency. In contrast to the HFB calculations, low-lying excited states in our 
approach are the RPA excitations (phonons) built on the vacuum states. Our vacuum states are yrast line states, i.e., 
the lowest energy states at a given rotational frequency. Note, that RPA phonons describe collective and noncollective 
excitations equally [3, E^l • The rotational bands are composed of the states with a common structure (characterized 
by the same parity, signature and connected by strong B(E2) transitions). We will demonstrate that the positive 
signature excitations are closely related to the shape transitions, that take place in the considered nuclei undergoing 
backbending. Hereafter, for the sake of discussion, we call our approach the GRPA. The validity of our approach will 
be confirmed by a remarkable agreement between available experimental data and our results for various quantities 
like kinematical and dynamical moment of inertia, quadrupole transitions etc. 

The paper is organized as follows: in Section H we review the Hartree-Bogoluibov approximation for rotating nuclei 
and discuss mean held results. Section HI is devoted to the discussion of positive signature RPA excitations and their 
relation to the backbending phenomenon. The conclusions are hnally drawn in Sec. IV. 
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II. THE MEAN FIELD SOLUTION 
A. The Hartree-Bogoliubov approximation 


We start with the Hamiltonian 

Hq = H - nnii = Ho- XrNr - nnh + v 

T—n,p 

= Ho- hflli + V (1) 

The unperturbed term consists of two pieces 
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The first is the Nilsson Hamiltonian 
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is a triaxial harmonic oscillator (HO) potential, whose frequencies satisfy the volume conserving condition ujiuj 2 td 3 = Wq 
{hujo = 41 MeV). In the cranking model with the standard Nilsson potential Q the value of the moment of 

inertia is largely overestimated due to the presence of the velocity dependent I ^ term. This term favours s.p. orbitals 
with large orbital momenta I and drives a nucleus to a rigid body rotation too fast. This shortcoming can be overcome 
by introducing the additional term. The second piece of Ho restores the local Galilean invariance broken in the rotating 
coordinate system and has the form 


hadd = XlmuJooK 


2 [r^Sx — xr- s) 
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The two-body potential has the following structure 


(5) 


V = Vpp + Vqq + VmM + Wcrcr- ( 6 ) 

It includes a monopole pairing, Vpp = —J2 t=p n^-^^T-Pr, where index y is labelling a com¬ 

plete set of the oscillator quantum numbers (|x) = \Nljm)) and the index y denotes the time-conjugated state [^. 
Vqq and Vmm are, respectively, separable quadrupole-quadrupole, Vqq = -^Y^T=o,i^('^)J2r=±J2tj.=o,i,2^Qi^[r])^y 
and monopole-monopole, Vmm = X]t=o i potentials. Vaa is a spin-spin interaction, V^a = 

— 4 1 Kcr(ir) X)r=± recall that the K quantum number (a projection of the angular mo¬ 

mentum on the quantization axis) is not conserved in rotating nonaxially deformed systems. However, the cranking 
Hamiltonian cj adheres to the D 2 spatial symmetry with respect to rotation by the angle tt around the rotational 
axis Xi- Consequently, all rotational states can be classified by the quantum number called signature r = exp(—wa) 
leading to selection rules for the total angular momentum I = a 2n, n = 0, ±1, ±2... In particular, in even-even 
nuclei the yrast band characterized by the positive signature quantum number r = -|-I (a = 0) consists of even spins 
only. All the one-body fields have good z component of isospin operator G and signature r. Multipole and spin- 
multipole fields of good signature are defined in Ref. . The tilde indicates that monopole and quadrupole fields are 
expressed in terms of doubly stretched coordinates Xi = [uji/ujo) Xi [^ . 

Using the generalized Bogoliubov transformation for quasiparticles (for example, for the positive signature quasi¬ 
particle we have + '^kVk) variational principle (see details in Refl^. we obtain the Hartree- 

Bogoliubov (HB) equations for the positive signature quasiparticle energies Ej (protons or neutrons) 


( h{l) A 
I A h{2) 


U, 



( 7 ) 
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Here, h{l)ki = iHo)ki - ^{h)kh {h{2))ki = -{Ho)ki - ^{Ii)ki and Aki = SkiGr < Pr > and \k > denotes a s.p. 
state of a Goodman spherical basis (see Ref|3l|). It is enough to solve the HB equations for the positive signature, 
since the negative signature eigenvalues and eigenvectors are obtained from the positive ones according to relation 

^ ( 8 ) 


where the state i denotes the signature partner of i . For a given value of the rotational frequency Q the quasiparticle 
(HB) vacuum state is defined as ai\) = aj|) = 0. 

The solution of a system nonlinear HB equations 0 is a nontrivial problem. In principle, the pairing gap should 
be determined self-consistently at each rotational frequency. However, in the vicinity of the backbending, the so¬ 
lution becomes highly unstable. In order to avoid unwanted singularities for certain values of H, we followed the 
phenomenological prescription 


Ar{n) 


A ,( 0 ) i (%)2 


(9) 


where He is the critical rotational frequency of the first band crossing. 

In general, in standard calculations with the Nilsson or Woods-Saxon potentials, the equilibrium deformations are 
determined with the aid of the Strutinsky procedure Q- The procedure, being a very effective tool for an analysis of 
experimental data related to the ground or yrast states, produces deformation parameters that are slightly different 
from those of the mean field calculations. The use of the former parameters (based on the Strutinsky procedure) 
for the RPA violates the self-consistency between the mean field and the RPA description. Therefore, to keep a 
self-consistency between the mean field and the RPA as much as possible, we use the recipe described below. 

It is well known mmiiiisi that for a deformed HO Hamiltonian, the quadrupole fields in double-stretched 
coordinates fulfill the stability conditions 

(Q^)=0, m = 0,1,2 (10) 

if nuclear self-consistency 

u;l{xl)=Lol{xl)=Lol{xl) ( 11 ) 


is satisfied in addition to the volume conserving constraint. Here, < ... > means the averaging over the mean field 
vacuum state of the rotating system. In virtue of the stability conditions the interaction will not distort further 
the deformed HO potential, if the latter is generated as a Hartree field. To this purpose, one starts with an isotropic 
HO potential of frequency wq and, then, generates the deformed part of the potential from the (unstretched) separable 
quadrupole-quadrupole (QQ) interaction. The outcome of this procedure is 

2 2 
TTtUJ^ 

Vho = —^-towo/3cos7Qo[+] - mulf]sinjQ2[l] (12) 

where one can use the following parameterization of the quadrupole deformation in terms of /3 and 7 (see, for example, 
Refi): 


towq/3cos7 = k[0](Qo[+]) 
muQ^sinj = -k[0 ](Q2[+])- 
The triaxial form given by Eq. 0 follows from defining 


/ 5 Stt 

1 - 2f3\ —cos{'j - —i) ,J = 1,2,3 

V 47r 3 


(13) 


(14) 


Here, we follow the convention on the sign of 7 -deformation accepted in Ref|2- The Hartree conditions have the form 
given by Eqs. OS only for a spherical HO potential plus the QQ forces. Quite often, Eqs. are considered as self- 
consistent conditions for pairing-|-QQ model interaction. In practice, the use of these conditions is based on AN = 0 
mixing in a small configuration space around the Fermi energy, which includes only 3 shells. This restriction limits 
a description of physical observables like a mean field radius and vibrational excitations. Furthermore, if AN — 2 
mixing is included, the RPA correlations are overestimated In addition, the QQ-forces, without the volume 
conservation condition, fail to yield a minimum for a mean field energy of rotating superdeformed nuclei |^. Due 
to all these facts, we allow small deviations from Eqs. m and enforce only the stability conditions m, which are 
our self-consistent conditions for the mean field calculations. These, in fact, hold also in the presence of pairing (see 
below) and ensure the separation of the pure rotational mode from the intrinsic excitations for a cranked harmonic 
oscillator 1 ^. 
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B. Some HB results 


As was mentioned in Introduction, for our calculations we choose and There are enough available 

experimental data on spectral characteristics and electromagnetic decay of high spins in these nuclei |35l| . It is also 
known that these nuclei possess axially symmetric ground states and exhibit the backbending behaviour at high spins. 
Moreover, nuclei with Z~ 66 and N~ 90 attract a theoretical attention for a long time, since the cranking model 
predicts that high-j quasiparticles drive rotating nuclei to triaxial shape . 



FIG. 1: Equilibrium deformations in jS-'y plane as a function of the angular momentum I = (/i) — 1/2 (in units of h). 

We perform the following calculations: 

1. The Hamiltonian o includes the term hadd 0 that restores the local Galilean invariance broken in the Nilsson 
potential. This calculation will be further refereed as the calculation I. 

2. The Hamiltonian o does not include the term hadd- Since this term is responsible for the correct behaviour of 
the moment of inertia, such calculations provide the answer about its importance, for example, for the value of 
the band crossing frequency. All microscopic results reported in literature (excluding the analysis of octupole 
excitations [^) do not include such term in the Nilsson potential. We will refer this calculation as the calculation 

II 

Parameters of the Nilsson potential were taken from Ref. 113 • These parameters were determined from a systematic 
analysis of the experimental s.p. levels of deformed nuclei of rare earth and actinide regions. In our calculations we 
include all shells m to A^ = 9 and this configuration space was sufficient for the fulfilling 99% of the sum rule for E2 
strength (see Refill). In contrast to the analysis of Refs ls^lifjldll based on ’’single stretched” coordinate method that 
involves the AAI = 2 mixing only approximative, we take into account the AA = 2 mixing exactly. This improves the 
accuracy of the mean field calculations. For the values of the pairing gaps At.( 0) at zero rotational frequency we use 
the results of ReflU AAr(O) = 0.940 MeV, Ap(0) = 0.985 MeV for i^soy and AAr(O) = 0.967 MeV, Ap(0) = 0.975 
MeV for ^®^Yb. These values have been obtained in order to reproduce nuclear masses. For quantization of angular 
momentum we use the equation < Ii >= 1 + 1/2. Here, the term 1/2 is due to the Nambu-Goldstone mode that 
appears in the RPA (see Hefl^. 

As shown in FiglU the triaxiality of the mean field sets in at the rotational frequency He which triggers a backbending 
in the considered nuclei. We obtain the critical rotational frequencies, at which the first band crossing occurs, h^c = 
0.250, 0.301 MeV for ^®^Yb, ^^®Dy, respectively. The parameters so determined yield results in a better agreement 
with experiments, compared to the ones obtained in Ref. |36| for N ^ 90. Moreover, our equilibrium deformations 
are short from being the self-consistent solutions of the HB equations. The doubly stretched quadrupole moments 
(Qo[+]) ctnd (Q 2 [+]) are approximately zero for all values of the equilibrium deformation parameters, consistently 
with the stability conditions m (see FigEJ. Indeed, any deviation from the equilibrium values of the deformation 
parameters /3 and 7 results into a higher HB energy. The ’’double stretched” monopole moment (Mo[°]) is not far 
from the corresponding standard one. The both values, the standard and the ’’double stretched” monopole moments, 
are almost independent on H. We infer from the just discussed tests that our solutions are close to the self-consistent 
HB ones. In contrast with our results, in Refl36l the fixed parameters (deformation, pairing gap) were used for all 
values of the rotational frequency. In addition, the analysis of Refl^ was based on the fitted moments of inertia, that 
were kept constant for all rotational frequencies. It is evident that such an analysis can be only used for a discussion 
of a general trend. 

We get more insight into the backbending mechanism if consider the potential energy surface near the transition 
point. As is shown in FigEl the potential energy surface of the total mean field energy i?o(/3, 7 ) = (Hq), for ^®®Dy at 
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FIG. 2: (Color online) The rotational behaviour of the calculated monopole and quadrupole moments. The ’’double stretched” 
and standard values are connected by dashed and solid line, respectively. 


hfl = 0.300 MeV (before the bifurcation point) exhibits the minimum for the axially symmetric shape which is lower 
than the minimum for strongly triaxial shape with 7 « 20°. The increase of the rotational frequency breaks the axial 
symmetry and a nucleus settles at the nonaxial minimum at hil = 0.302 MeV (after the bifurcation point). Notice, 
the main difference in these minima is a strong nonaxial deformation of the one minimum in contrast to the other, 
while the /3- deformation is almost the same. 

Dealing with transitional nuclei, however, the minimum becomes very shallow for a collective (around xi rotational 
axis) and non-collective (around xy, symmetry axis) rotation as the rotational frequency increases. In fact, the energy 
minima for the collective (« —1730.6 MeV) and non-collective (« —1730.62 MeV) rotations are almost degenerate 
near the crossing point of the ground with the 7 - band for ^®°Dy (see the upper panel of the Fig. EJ). The energy 
difference is about 15 keV near the critical rotational frequency where the backbending occurs. At the transition 
point, the competition between collective and non-collective rotations breaks the axial symmetry and yields nonaxial 
shapes. Does this behaviour correspond to a phase transition ? 

First, notice that a half of an experimental value for 7 -vibrations hD,K=2/‘2, = 0.414 MeV (at HD, = 0) is 
close to the collective rotational frequency hilc = 0.301 MeV at which the shape transition occurs. Second, let us 
consider an axially deformed system, defined by the Hamiltonian H in the laboratory frame, that rotates about a 
symmetry axis z with a rotational frequency D. The angular momentum is a good quantum number and, consequently, 
[Jz,0\^] = ATOjf. Here, the phonon describes the vibrational state with K being the value of the angular 
momentum carried by the phonons along the symmetry axis, z axis. Thus, one obtains 

[Hn,Ol] = [H- = {WK - = lokO^j,, (15) 

where ujk is the phonon energy of the mode K in the laboratory frame at D = 0. This equation implies that at the 
rotational frequency D,cr = ujk/K one of the RPA frequency uk vanishes in the rotating frame (see discussion in 
Refsj^ EJi Ed)- At this point of bifurcation we could expect the SSB effect of the rotating mean field due to the 
appearance of the Goldstone boson related to the multipole-multipole forces with quantum number K. For an axially 
deformed system, one obtains the breaking of the axial symmetry, since the lowest critical frequency corresponds to 
7 -vibrations with K = 2 mil. The rotation around collective xi axis couples, however, vibrational modes with 
different K and the critical rotational frequency is influenced by this coupling: it becomes lower. The value of D,cr can 
be affected by the degree of the collectivity of the vibrational excitations, as we will see below. The most important 
outcome from this consideration is that in the vicinity of the shape transition there is an anomalously low vibrational 
mode related to the deformation parameter 7 . 

For the shape transition takes place at the rotational frequency TiDc « 0.25 MeV, while the experimental 

bifurcation point (a half of the 7 -vibrational excitation energy at Tiff = 0) is TiDcr ~ 0.45 MeV. The energy difference 
at the transition point between the collective (~ —1837.7 MeV) and the non-collective (~ —1836.6 MeV) minima 
is still large ^ 1 MeV (see the lower panel of FigEJ. As one can see in FigEl the difference between the axially 
symmetric minimum at hD. = 0.245 MeV (before the transition) and the nonaxial one at = 0.255 MeV (after 
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FIG. 3: (Color online) The potential energy snrface Eq{P, 7 ) for before (TlQ, = 0.300 MeV) and after {hQ = 0.302 MeV) 

the transition point. There are two local minima with almost the same 13 deformation and completely different 7 -deformation. 



tiQ [MeV] 

FIG. 4: (Color online) The rotational dependence of the total mean field energy i?n(/3,7 = 0) = (Hn) for the axially symmetric 
equilibrium deformation. Results for collective {xi axis) and noncollective (*3 axis) rotation are connected by solid and dashed 
line, respectively. Results for ^®®Dy and '^®^Yb are displayed on upper and lower panel, respectively. 


the transition) for a collective rotation is about ~ 1.2 MeV. The deformation parameters change smoothly with the 
increase of the rotational frequency at the transition point (see FigOJ. It appears, that for there is a different 

mechanism responsible for the observed backbending. 

To elucidate the different character of the shape transition from axially symmetric to the triaxial shape and its 
relation to a phase transition, we consider potential landscape sections in the vicinity of the shape transition. The 
phase transition is detected by means of the order parameter as a function of a control parameter 0. In our case, 
the deformation parameters (3 and 7 are natural order parameters, while the rotational frequency 17 is a control 










































FIG. 5: (Color online) The potential energy surface En{l^,'y) for before {HQ = 0.245 MeV) and after (hO, = 0.255 MeV) 

the transition point. 


parameter that characterizes a rotational state in the rotating frame. Since we analyze a shape transition from the 
axially symmetric shape (7 = 0 ) to the triaxial one (7 yf 0 ), we choose only the deformation parameter 7 as the order 
parameter that reflects the broken axial symmetry. Such a choice is well justified, since the deformation parameter (3 
preserves its value before and after the shape transition in both nuclei: f3t « 0.2 for and f3t « 0.31 for 

Thus, we consider a mean field value of the cranking Hamiltonian, f3t) = {Hq), for different values of H (our 

state variable) and 7 (order parameter) at fixed value of Pt ■ 

For we observe the emergence of the order parameter 7 above the critical value Tiflc = 0.301 MeV of the 

control parameter Q (see a top panel in Fig0. Below and above the transition point there is a unique phase whose 
properties are continuously connected to one of the coexistent phases at the transition point. The order parameter 
changes discontinuously as the nucleus passes through the critical point from axially symmetric shape to the triaxial 
one. The polynomial fit of the potential landscape section at fi^lc = 0.301 MeV yields the following expression 


i^(H;7) = Fo(fl) + F2{n)-/^ - + F^n)^^ 


where the coefficients Fo{Q) = 0.3169 MeV, 7 in degrees and F 2 {il.) = 0.12239, ^ 3 ( 0 ) 
are defined in corresponding units. We can transform this polynomial to the form 


where 


F = 


F(H;7)-fo(») 

Fn 


rf rf 

“y ~ y 


iF_ 

4 


Fn = 


(3^3)^ 

{mf 


8F2F4 


4 F 4 


(16) 

0.009199, F 4 (H) = 1.7 X 10-4 

(17) 


(18) 


The expression represents the generic form of the anharmonic model of the structural first order phase transitions 
in condensed matter physics (cf [dg). The condition dF/drj = 0 determines the following solutions for the order 
parameter rj 


1 ± VI - 4a 

T] = 


7 = 0, 


2 


(19) 
















9 


156 


50 

40 

30 

20 

10 

0 


Dy 


0 

-10 

-20 

-30 

-40 

-50 


' I ' I ' I I I ' 

0 0.2 0.4 

m[Mey] 

162 

Yb 


0 0.2 0.4 

m[Mey] 


tiQ = 0.300 MeV 

3.0 

? 2.0 

— 1.0 

U. 

0.0 


0.301 MeV 0.302 MeV 




I I I I I I I I I I I I I 

0 10 20 




tiQ = 0.245 MeV 

1.2 


0.250 MeV 0.255 MeV 


0.8 

0.4 

0.0 



-8 -4 0 4 8 

Y 


FIG. 6 : (Color online) The rotational dependence of the order parameter 7 and the energy surfaces sections F{Q,,'y) = 
E{'y,l3t) — Emin for Dy (top) and (bottom) before and after the transition point. The energy is given relative to the 

value Emin = EniPt,^) at hfl = 0.255, 0.302 MeV for respectively. 


If a > 1/4, the functional F has a single minimum at 77 = 0. Depending on values of a, defined in the interval 
0 < a < 1/4, the functional F manifests the transition from one stable minimum at zero order parameter via 
one minimum+metastable state to the other stable minimum with the nonzero order parameter. In particular, at 
the universal value of a = 2/9 the functional F has two minimum values with F = 0 atr 7 = 0—and 
77 = 2/3 ^ 7 « 27*^ and a maximum at 77 = 1/3 ^ 7 « 13.5°. The correspondence between the actual value 7 ~ 20° 
and the one obtained from the generic model is quite good. Thus the backbending in possesses typical features 

of the first order phase transition. 

In the case of the energy and the order parameter (FigEll are smooth functions in the vicinity of 

the transition point Qc- This implies that two phases, 7 = 0 and 7 7^ 0, on either side of the transition point should 
coincide. Therefore, for Q near the transition point f 2 c we can expand our functional F{n, 7) = f3t) — Emin (see 

FigEJ in the form 


F{n; 7 ) = Fi(r !)7 + F 2 {n)-f^ + ^ 3 ( 11 ) 7 ° + + ... 

The conditions of the phase equilibrium (further we restrict the expansion 1201 ) up to the terms with 7 ^) 

dF 


97 


= Fi(fl) + F 2 (D )27 + F3(D)37^ + F 4 (11)47'' = 0 


d'^F 


= 2F2(D) + 6 F 3 (D )7 + 12^4(11)72 > 0 , 


which should be valid for all values of fl and 7 (including 7 = 0 ), yield 

Fi{n) = 0 

Eqs. 1211 ) and 1221 ) can be rewritten as 

2 F 2 (fI )7 + 3^3(11)72 + 4 E 4 (D) 7 ° = 0 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 


2 ^ 2 ( 0 ) + 6 ^ 3 ( 11)7 + 12 ^ 4 ( 11)72 > 0 


( 25 ) 
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which implies the following inequality 


2F4{n)-/^ > F2{n) 


(26) 


This inequality holds for all values of 7 (including 7 = 0 at = flc), which leads to F2{i2 = flc) < 0. On the other 
hand, from the stability condition Ea. (l 22 ll at the transition point Oc and 7 = 0, we also have ^ 2(0 = Oc) > 0. The 
both inequalities can coincide only when ^ 2(0 = Oc) = 0. Using the result Fi(r2c) = F 2 {flc) = 0 and the fact that all 
phases at the transition point should coincide, we obtain from Eq. eu that ^ 3(0 = ^lc) = 0. Assuming that ^3 = 0 
for all O, the minimum condition Eq. (|23 yields the following solution for the order parameter 


71 = 0, 


F2{n) 

2F4(0) 


^ 0 for O ^ Oc 
= 0 for Q = Qc 


(27) 


Since at the transition point F 2 (Uc) = 0, one can propose the following definition of the function ^ 2 ( 0 ): 

F2{n)^^^P- (n-n,) (28) 

ail 

Thus, we have 7 ~ (U — Uc)"^ and the critical exponent ly = 1/2, in accord with the classical Landau theory, where 
the temperature is replaced by the rotational frequency. 

Our extension of the Landau-type approach for rotating nuclei is nicely confirmed by the numerical results. The 
polynomial fit of the energy potential surfaces for (Fig0 yields Fi(0) = ^ 3 ( 0 ) = 0 for all considered values of 

the rotational frequencies and F 2 (Uc) = 0 at hUc = 0.25 MeV. Moreover, in the vicinity of flc we obtain dF2{fl)/hdFl k, 
—3.5 (^ 4 ( 0 ) > 0 for all U). In an agreement with Eos. 12711 . 12811 . we have only the phase 7 = 0 for TiFl < Mlc and the 
phase 7^0 for TiU, > 7if2c. The energy surfaces are symmetric with regard of the sign of 7 and this also supports 
the idea that the effective energy F can be expressed as an analytic function of the order parameter 7 . Thus, the 
backbending in can be classified as the phase transition of the second order. 


C. Quasiparticle spectra 

To understand the microscopic origin of the quantum shape-phase transitions, we analyse first the quasiparticle 
spectra and the rotational evolution of different observables like quadrupole moments and moments of inertia. A 
numerical analysis of the expectation value of the nonaxial quadrupole moment 

(Q 2 [^]) = E (WlW) E (29) 

kl i 

shows that in ^®^Yb high-j neutron and proton orbitals that belong to ii 3/2 and / 111/2 subshells, respectively, give 
the main contribution to the the expectation value. The nonaxial deformation grows due to the rotational alignment 
(RAL). The crossing frequencies, where the configurations with aligned quasiparticles become yrast, are order of 
Qc ~ 2A/i, where i is an aligned angular momentum carried by quasiparticles. Since the neutron gap is smaller than 
the proton gap, one may expect that the backbending should occur due to the alignment of the quasineutron orbitals 
that could contribute i ~ 8h. 

We trace the rotational evolution of quasiparticle orbitals in the rotating frame (routhians) as a function of the 
equilibrium parameters (e, 7 , A). At Q = 0 each orbital is characterized by the asymptotic Nilsson quantum numbers. 
However, these numbers lose their validity in the rotating case due to a strong mixing. Hereafter, they are used only 
for convenience. The analysis of the routhians for neutrons (Fig[7I) and for protons (FigO indicates that the lowest 
quasicrossings occur: at hil 0.275(0.42) MeV for neutron (proton) system in ^^®Dy; at hQ 0.245(0.41) MeV 
for neutron (proton) system in ^®^Yb. We recall that the shape-phase transition occurs at Ti^lc ~ 0.25,0.3 MeV in 
i62Yb^ i56Dy, respectively. The proximity of the critical point to the two-quasiparticle neutron quasicrossing in both 
investigated nuclei (especially, in ^^^Yb) suggests that the alignment of a pair ^ 13/2 is the main mechanism that drives 
both nuclei to triaxial shapes (cf pM This mechanism itself, however, does not provide the explanation for the 
different character of the shape-phase transition. We recall that the routhians exhibit the dynamics of noninteracting 
quasiparticles. Evidently, the interaction between quasiparticle orbitals is important and could change substantially 
as a function of the neutron and proton numbers. Indeed, as we will see below (see Sec.HI), this explains the type of 
shape-phase transitions discussed in the present section. 
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FIG. 7: (Color online) Lowest quasineutron energies for (upper panel) and (lower panel). Thick (thin) lines are 

used for the positive (negative) parity states. The positive (negative) signature states are connected by solid (dashed) lines. 
At n = 0 the levels A, B,C, D, E correspond to the Nilsson states: 3/2[651] (the subshell ii 3/2 ), l/2[660]) (the subshell 113 / 2 ), 
3/2[521] (the subshell / 19 / 2 ), 5/2[521](the subshell fr/ 2 ), ll/2[505] (the subshell / 111 / 2 ), respectively. The shape transition point 
is denoted by the vertical line. The quasicrossing is surrounded by a circle. 


D. Inertial properties 

The moment of inertia is the benchmark for microscopic models of nuclear rotation. To understand the interplay 
between s.p. degrees of freedom and collective effects we calculate the kinematical = {ii)/n and dynamical 

/o\ ^ 

3) ' = d{Ii)/d£t moments of inertia and compare with the corresponding experimental values. The kinematical 
moment of inertia reflects the collective properties of the rotating mean field. The dynamical moment of inertia 
due to obvious relation /dfl is very sensitive to structural changes of the mean field. It reflects the 

rearrangement of the two-body interaction upon rotation which leads to level crossings and shift in the deformation. 
In fact, the dynamical moment of inertia provides a definite criteria on the self-consistency of the rotating mean field 
when it is compared with the Thouless-Valatin moment of inertia calculated in the RPA. A full self-consistency is 
achieved if they are equal (see discussion in Ref. |^). 

In Figl^the experimental and theoretical values of kinematical (upper panels) and dynamical (lower panels) mo¬ 
ments of inertia are compared for ^®^Yb (left panels) and ^^®Dy (right panels). We remind that the calculations 
I (II) include (not include) the term hadd, Eq.@. While both calculations reproduce the rotational evolution of 
the kinematical moment of inertia, the agreement with the experimental data is much better for calculations I. It 
is interesting to note that with the increase of the rotational frequency the ratio increases from ~ 50% at 

hil ^ 0.2 MeV to ^ 85% at hil ^ 0.5 MeV. Partially, the effect of broken Galilean symmetry is reduced due to the 
alignment of the high-j intruder states with a large orbital momentum 1. These states contribute to the collective 
angular momentum and decrease the effect of the I ^ term in the Nilsson potential. 

Although one observes a similar pattern for the backbending in the considered nuclei (upper panels, FigEI), a different 
response of a nuclear field upon the rotation becomes more evident with the aid of the experimental dynamical moment 
of inertia = dl/dfl ~ A/AE^ as a function of the angular frequency (see the second panel from the bottom, 
Fig0. Here, hil = Ejl2, Ej is the 7 -transition energy between two neighboring states that differ on two units of the 
angular momentum / and AE-^ is the difference between two consecutive 7 -transitions. At the transition point 
wildly fluctuates with a huge amplitude in whereas these fluctuations are quite mild in This behaviour 
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FIG. 8: (Color online) Lowest quasiproton energies for (upper panel) and (lower panel). Thick (thin) lines are used 

for the negative (positive) parity states. At 12 = 0 the levels A,B,C,D correspond to the Nilsson states: 7/2[523] (the subshell 
ftii/ 2 ), 5/2[532] (the subshell hii/ 2 ), 3/2[411](the subshell d^/ 2 ), 5/2[413] (the subshell 57 / 2 ), respectively, in ^®®Dy; 7/2[523] 
(the subshell hii/ 2 ), 9/2[514] (the subshell hii/ 2 ), 5/2[402] (the subshell ^ 5 / 2 ), 7/2[404] (the subshell 37 / 2 ), respectively, in 
i62Yb 'pjje shape transition point is denoted by the vertical line. The quasicrossing is surrounded by a circle. 


can be understood by virtue of the expansion of the microscopic moment of inertia at small rotational frequency 

J7b+Jm + ^^|§ + ... (30) 

The Inglis-Belyaev moment of inertia, fT/sj neglects a residual two-body interaction between quasiparticles. Its 
behaviour is similar to the one of the kinematical moment of inertia, . The second term JTm is a Migdal moment 
of inertia |49|. resulting from the effect of rotation on the residual two-body interaction and, in particular, on the 
pairing interaction. In our calculations the pairing gaps change smoothly in accord with the phenomenological 
prescription, Eq. ©■ The third term describes the variation of the self-consistent mean field, namely, the change of 
the deformation (/3i.2 = /?, 7 ) under rotation. Therefore, a drastic change of the mean field configuration (namely, 
7 -deformation) in ^^®Dy (see Fig|^ top panels) explains large fluctuations of the dynamical moment of inertia at 
the transition point. In contrast, the smooth behaviour of the function F at the transition point (see FigEl bottom 
panels) implies a small amplitude of fluctuations of the dynamical moment of inertia in ^®^Yb. The magnitude of 
fluctuations can be traced by means of the ratio 11 = as a function of the rotational frequency (see the 

bottom panel in FiglHI). While this ratio is about 11 ^ 2 at the transition point in ^®^Yb, in ^®®Dy it is much larger 
n ~ 8 . Referring to above analysis of the mean field solutions, we can formulate the empirical rule to detect the order 
of the quantum shape-phase transition at the backbending: if the ratio 11 = >> 1 at the transition point, 

the shape transition can be associated with a first order phase transition. 

The nature of the backbending becomes more evident by dint of the RPA analysis presented below. For a com¬ 
pleteness, we also compare the Thouless-Valatin moment of inertia (the definition is presented in Sec.Ill) with the 
dynamical moment of inertia calculated in the mean field approximation. Notice, that all three terms in Eq. 
are included into the Thouless-Valatin moment of inertia, ^rpa = O', calculated in the RPA (see, for example, the 
discussion in Pefl^ for the exact model without pairing), which is valid at large rotational frequencies as well. One 
can observe a remarkable agreement between experimental, mean field and CRPA results. The agreement between 
the mean field and CRPA results confirms a good accuracy of our mean field calculations. This agreement, on the 
other hand, demonstrates a validity of our CRPA approach in the backbending region, at least, for the description of 











13 



FIG. 9: (Color online) The kinematical (top panels) and dynamical Q'P^n) (second panel from the bottom) moments 

of inertia (in units Tl^/M eV) are compared with the corresponding experimental values (filled squares). Experimental values are 
connected by dash-dotted line to guide eyes. The results for the kinematical moment of inertia Q'j^^(fl), with (the calculations 
I) and without (the calculations II) the additional term hadd, Eq. 01 are connected by solid and dashed lines, respectively. 
The second panel (from the top) displays the ratio R = for each nucleus. The experimental and calculated dynamical 

moment of inertia (solid line) (obtained by the calculation I) are compared with the Thouless-Valatin inertia moment 
Qrpa , Eg. 18911 (dashed line). In the bottom panels the ratio, 11 = is displayed. 


the spectrum. 


III. QUADRUPOLE COLLECTIVE EXCITATIONS 
A. Quasiparticle RPA in rotating systems 

In this section we discuss the RPA results to get a further insight into the backbending mechanism. By means 
of a generalized Bogoliubov transformation, we express the Hamiltonian given by Eq. o in terms of quasiparticle 
creation {a\) and annihilation (a^) operators. We then face with the RPA equations of motion, written in the form 

mim 

[Ha,Pi,] = ihujlX^, [Ha,X^] = -ihP„, 

[X,,P,,] = ihS,,,, (31) 

where Xi,, Pi, are, respectively, the collective coordinates and their conjugate momenta. The solution of the above 
equations yields the RPA eigenvalues huj„ and eigenfunctions 

|i.>= Oi\RPA> = [^X„ - 

v2 ^ V “ ' 


( 32 ) 
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?.+ _ ^.+^,+ ht, = -+ 


ttil at, b^j = a^a^ are used. The first 


where fi = (kl) or (kl, kl). Here, the boson-like operators — u^. 

equality introduces the positive signature boson, while the other two determine the negative signature ones. These 
two-quasiparticle operators are treated in the quasi-boson approximation (QBA) as an elementary bosons, i.e., all 
commutators between them are approximated by their expectation values with the uncorrelated HB vacuum Q. The 
corresponding commutation relations can be found in Refl25l In this approximation any single-particle operator F 
can be expressed as F = (F) + F^^'> + where the second and third terms are linear and bilinear order terms in 
the boson expansion. We recall that in the QBA one includes all second order terms into the boson Hamiltonian such 
that {F — (F))^ = The positive and negative signature boson spaces are not mixed, since the corresponding 

operators commute and F[q = F{fi{r = -|-1) -|- Flair = —1). Consequently, we can solve the eigenvalue equations Id II) 
for Ha{+) and Fo(—), separately. 

The symmetry properties of the cranking Hamiltonian yield 


[Hal+), NT^n,p]RPA — 0 , [Hal+), — 0. 

The presence of the cranking term in Hamiltonian o leads to the following equation 

[Hnl-),rt = nnr\ 


where Pf = (/2 + ih )/and T = (ri)i = (J 2 — il ^)/fulfill the commutation relation 

[r,rt] = A 


(33) 


(34) 


(35) 


According to Eqs. P|) . we have two Nambu-Goldstone modes, one is associated with the violation of the conservation 
law for a particle number, the other is a positive signature zero frequency rotational solution associated with the 
breaking of spherical symmetry. Eq. El, on the other hand, yields a negative signature redundant solution of energy 
Wi, = H, which describes a collective rotational mode arising from the symmetries broken by the external rotational 
field (the cranking term). Eqs. (TTHll and (IHH) ensure the separation of the spurious or redundant solutions from 
the intrinsic ones. They would be automatically satisfied if the single-particle basis was generated by means of 
a self-consistent HB calculation. As we shall show, they are fulfilled with a good accuracy also in our, not fully 
self-consistent, HB treatment. 

We recall that the yrast states possess the positive signature quantum number. Obviously, SSB effects of the mean 
field are related to Nambu-Goldstone modes of the same signature. Therefore, in this paper our analysis is focused 
upon positive signature RPA excitations. 

The positive signature Hamiltonian consists of the following terms 




r=iV,P 
1 


- 1/2 ^ Y. 


T=0,1 


T=0,1 


- 9 E « 2 [T] ^ 


(36) 


T=0,1 


M=0.1,2 


Here, = Si -\- ej are two-quasiparticle energies and the definition of matrix elements of the operators involved in 
the Hamiltonian II36I) can be found in Ref l^ The solution of the equations of motion, Eq. EU, leads to the following 
determinant of the secular equations 


Flu}^) = det (R- — ), 
with a dimension n = 12 and c = ko,K 2 or Gt- The matrix elements Rkml^v) = X) 


(37) 

Ct^/lEl - ojD 


fi QkpQmp^fj, 

involve the coefficients or 01 ^, for different combinations of matrix elements qkp (see details in Refs l25ll31ll 

The zeros of the function Tlujy) = 0 yield the CRPA eigenfrequencies Once the RPA solutions are found, the 
Hamiltonian EEl can be written in terms of the collective modes (A,^, P^) 

1 
2 








( 38 ) 
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where gi^ = l/’^SB.PA^ gN^ are the Thouless-Valatin moment of inertia and the nucleus mass, respectively, determined 
in the uniform rotating (UR) frame. The general derivation of the mass parameters gj-^ and g^^ can be found in 
RefsIHIE^ In particular, for the moment of inertia we have 


'^RPA — 


detA 

detB 


(39) 


where A^i is the matrix (the dimension n = 10 ) given by the part of the matrix !Fki{ui = 0 ) corresponding to the s.p. 
operators k, and I involved in H[r = +1] with the structure ( 6 it type. The matrix B^i (with the dimension 

fWj(i) 

n = 11 ) is given by matrix Aki with additional one column and row involving 0 and sums ± e 

where are quasiparticle matrix elements of the operator Ji (see Reflhol for details). 

Transition probabilities for XX transition lip \I'v' > between two high-spin states is given by expression 


B{X\]Iv ^ I'v') = (40) 

{II\gi\I' I'f\ < p'\M^^\XX-gi = I' - I)\v > p 

At high spin limit (/ >> A, /' >> A), this expression takes the form 


B{XX-Iv I'yr) ~ 

I < RPA\ \Q^,M^^\XX-gi= I' - I)] \RPA> p 


(41) 


for the transition from one phonon states into the yrast line states. Here, {XXgi) is the linear boson part of the 

corresponding transition operator of type X, multipolarity A and the projection gi onto the rotation axis (the first axis 
of the UR system). The commutator in (ITO can be easily expressed in terms of phonon amplitudes (p = U 

or and | RPA) denotes the RPA vacuum (yrast state) at the rotational frequency U. The multipole operator in 
the rotating frame is obtained from the corresponding one in the laboratory frame according to the prescription M 

MiXXy,)=J2 A4(AAm3). (42) 

M3 


Taking into account that < J^|.^ 2 m 3 =o, 2 [+]I^ >= (l-^ 2 ^ 3 =o, 2 [+]l) holds in the first RPA order and Ea. (l^ . we have 


r(B) 


B{E2;Ip ^ I -2p) = < {E2-, m = 2)\v > 






(E) 


(43) 


Using the /3 - 7 deformation parameterization, Eq. dia, and the oscillator value for the isoscalar quadrupole constant, 
K 2 [ 0 ] « ^ ~ ^mu^R? {R = 1.2A^/^/to), we obtain an approximate expression 


B{E2;Iyr^ I—2yr)pi 

1 4 2 2t'E N 

cos^(- - 7 ) 


2 167r2 


'6 


(44) 


Here, we use = {eZ/A)M. From this expression one can conclude that the change of the 7 -deformation from 

zero to 7 r /6 will increase the transition probability. The negative sign of 7 -deformation results in the decrease of this 
probability. As we will see below, this expression is useful for the analysis of the experimental data. 


B. Determination of the residual interaction strengths 


In order to determine the strength n of the monopole and quadrupole interactions, we start with the oscillator 
values 


«a[0] 


47r TOWg 

2A-blA<f2A-2 


«a[ 1 ] 


ttVi 


(45) 


A < f2'’' > 
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FIG. 10: (Color online) The ratio / = Kcai/K'osc between actual and oscillator isoscalar quadrupole strength constants for the 
positive and negative signature is displayed on the upper panel. On the lower panel, the actual and the BCS pairing constants 
Gjv and Gp (see the definition in the text) are shown. All actual values used in the calculations are obtained as to fulfill the 
RPA equations, Eqs. (OSJ-lEl, for the spurious or redundant modes. The results correspond to the calculation I, when the 
term hadd 0, that restores the local Galilean invariance broken in the Nilsson potential, is included. 


For instance, the isoscalar strength follows from enforcing the Hartree self-consistent conditions. We then change 
slightly these strengths and the pairing interaction constants Gt at each rotational frequency, while keeping constant 
the k[1]//c[ 0] ratio (ko[1]/ko[0] ~ —18, K 2 [ 1 ]/«: 2 [ 0 ] « —3.6), so as to fulfill the RPA equations (|^-l|33 for the spurious 
or redundant modes. 

The rotational dependence of these parameters is relatively weak. The ratio between actual and oscillator value 
for the isoscalar quadrupole constant is displayed in FiglTUl For comparison, the BCS values Gt, Gn ~ 22.blA 
MeV, Gp « 26.5/A MeV, obtained from the systematic of pairing gaps (see R ef l2^ . are shown in Fig^|as well. 
Surprisingly, the difference between actual values and the latter ones is mild. However, the determination of «; 2 [ 0 ] 
and Gt is a tedious task, since a tiny change of the strength parameters leads to large shifts in energies of spurious 
modes. The constants so determined differ from the HO ones by 5-10% at most. For the spin-spin interaction, we use 
the generally accepted strengths 

47r 

K^[0] = K„[l] = -28 —MeV 

J\. 

for all rotational frequencies. The spin-spin interaction does not influence the position of Nambu-Goldstone modes and, 
therefore, does not play any role in the self-consistent determination of the quadrupole strength constants. Finally, 
we adopted bare charges to compute the E2 strengths and a quenching factor gs = 0.7 for the spin gyromagnetic 
ratios to compute the Ml strengths. 

By using the above set of parameters, it was possible not only to separate the spurious and rotational solutions from 
the intrinsic modes, but also to reproduce the experimental dependence of the lowest (3- and 7 -bands on the rotational 
frequency hfl and, in particular, to observe the crossing of the 7 - band with the ground band in correspondence with 
the onset of triaxiality. 

Concluding the discussion about our calculation scheme, we should mention some drawbacks of our approach. One 
of the disadvantage is the CRPA breaks down at the transition point when Ap or A„ vanish [^. We have avoided 
this problem by means of the phenomenological prescription for the rotational dependence of the pairing gap (see 
Section H). In principle, projection methods should be used in the transition region in order to calculate transition 
matrix elements. In the case of the phase transition of the first order we run into the problem of the formidable 
overlap integral. A theory of large amplitude motion would provide a superior means to solve this problem. 

However, in contrast to the standard RPA calculations, where the residual strength constants are fixed for all values 
of ft (see e.g. [1323): determine the strength constants for each value of ft by the requirement of the validity of 

the conservation laws. This enables us to overcome the instability of the RPA calculations at the transition region, at 
least, for the excitations (see also discussion for quantum dots in R ef l2lll . Although the amplitudes (see Eq. (P|l) 
are higher for the RPA modes in the transition region than in other regions, the relation \ < \ is still valid to 



















17 


hold the QBA. And least but not last, the CRPA becomes quite effective at high spins, when the pairing correlations 
are suppressed. 


C. Positive signature excitations in the rotating frame 

In this section we analyse the spectral and decay properties of the positive signature states. Transition probabilities 
represent the most challenging task, since they are sensitive about the wave function structure of the considered state 
and, therefore, exhibit hidden drawbacks of a theory. 



FIG. 11: (Color online) Reduced transition probabilities B{E2\ I yr —> / — 2yr) along the yrast line. Experimental data 
(filled squares) are connected with a thin line to guide eyes. The results of calculations by means of En. M.SII and of Eq. are 
connected by solid and dashed lines, respectively. 

Experimental values of B{E2, I v ^ I' yr) are deduced from the half life of the yrast states using the standard, 
long wave limit expressions B{E2,i ^ /) = P{i /)/(!.223 x 10®A®) (e® fm"*) Here, the transition energy 

E^ is in MeV and the absolute transition probability P{i f) = In2/T(z ^ /) is the related to the half life 
T{i —> /) (in seconds). For ^^®Dy the yrast energies and corresponding half lifes are known up to the momentum 
I « 407i. For ^®®Yb yrast energies are observed up to / « 28fi but half lifes were measured only up to / « 207i. 
While the cranking approach should be complemented with a projection technique in the backbending region due 
large fluctuation of the angular momentum (cf RefH, its validity becomes much better at high spins. The agreement 
between calculated and experimental values of intraband B{E2) transitions along the yrast line is especially good 
after the transition point. 

Experimental data are compared with the results of calculations (a) by means of Eg. (I43II and calculations (b) by 
means of Eq. lEll (see Figtmi. In the calculations (a) we use the mean field values for the quadrupole operators. In the 
calculations (b) these values have been replaced by the deformation parameters via Eq. (ICT and we use the oscillator 
value for the quadrupole strength constant. The calculations (a) evidently manifest the backbending effect obtained 
for the moments of inertia (see Fig0 at hUc « 0.25, 0.3 MeV for ^®®Yb and ^®®Dy, respectively. Thus, the use of the 
self-consistent expectation values is crucial to reproduce the experimental behaviour of the yrast band 

decay. The calculations (b) (Eg. (Ell ) reproduce the experimental data with less accuracy. However, these results also 
catch on the correlation between the sign of the y-deformation and the behaviour of the transition probability. The 
onset of the positive (negative) values of y-deformation lead to the increase (decrease) of the transition probability 
along the yrast line. This fact nicely correlates with the experimental data. 

To analyse experimental data on low lying excited states near the yrast line we construct the Routhian function 
for each rotational band v {v = yrast, ...) 


R, = E,(H) - nniin), n{i) = + (46) 

and define the experimental excitation energy in the rotating frame = i?i,(H) — Ryr{t2) [^. This energy can be 

directly compared with the corresponding solutions huji, of the RPA secular equations for a given rotational frequency 


















18 


156 


V=1 


Dy 


v = 2 




hQ [ MeV ] 


FIG. 12: (Color online) Experimental (filed square) and calculated excitation energies in the rotating frame = 

Ri^{Q.) — Ryr{^) (upper panel) and reduced transition probabilities B{E2\Iv —> I'yr) (lower panel) for two lowest RPA 
solutions, V = 1 (left) and v — 2 (right), as a function of the angular frequency hSl. The results for excitations, calculated 
with and without the term h^dd, Eq. are connected by solid and dash-dotted lines, respectively. The experimental values 
for the transitions are indicated by arrows. The tansitions with AI = 2, 0, —2 are connected by solid, dash-dotted and dashed 
lines, respectively. Two-quasiparticle energies, indicated by arrows with numbers 1, 2, and 3 in circles, originate from two s.p. 
Nilsson states (nh |[651] |[651]), (pp|[523] |[523]), and (nh i[530] |[523]) at If = 0 and At = 0. 


We remind that our vacuum states are adiabatic quasiparticle configurations that correspond to the occupied 
orbitals at a given rotational frequency. 

The first lowest positive signature RPA solutions {v = 1) create an excited band built on the yrast line. In both 
nuclei, this band is close to the experimental y-excitations with even spin at small rotational frequency Q, (see FigsCU 
n^i . In the low-spin region, the transitions B{E2; Ji/ = 1 —> I'yr) are large and RPA solutions exhibit a strong 
collective nature of the 7 - band (Weisskopf units are B{E2)w = 49.9, 52.5 for respectively) in 

both nuclei. The collective character of the low-spin part of the 7 - bands is manifested in the phonon structure that 
is composed by a few two-quasiparticle components (see below and Tables I,II). 

In ^^®Dy the RPA results, obtained with the term hadd ®, reproduce quite well the rotational behaviour of the 
lowest excitations associated with the collective 7 -vibrations. We also obtain a good agreement with the experimental 
value of the critical rotational frequency hi^cr ~ 0.324 MeV, at which this mode disappear in the rotating frame. This 
results is very close to the critical rotational frequency hilc « 0.3 MeV (see FigEl, where the backbending occurs (see 
Ea. (ll5ll and discussion in Sec. II). In accordance with the phenomenological theory of the first order phase transitions 
discussed in Section II, this soft mode creates a shape transition from the axial to the nonaxial shapes. We emphasize 
that without the term hadd the transition point is located considerably higher than the experimental one. In ^®^Yb, 
the results with and without this term are very close (see Fig[T3|). In contrast to ^®®Dy, in ^®^Yb the crossing, is 
determined by a single two-quasiparticle state (see Table I at / = lOh). The term hadd results in a collective effect to 
the mean field solution, while its contribution to a single, two-quasiparticle energy is weak. As discussed in Sec. 
II, the effect of this term is also reduced due to the alignment. 

In Tables in miwe present the contribution = 1) = ~ main two-quasiparticle components 

to the norm = 1) = 1 of the RPA solution as a function of the angular momentum (rotational 

frequency). The collectivity is weaker in ^®^Yb than in ^®®Dy. With the increase of the rotational frequency, in ^®^Yb 
the phonon loses the collective nature and a single, two-quasiparticle neutron component {ij = nh|[651] |[651]) is 
dominant in the transition region. At the crossing point its weight reaches 96% and B{E2-, Iv = 1 ^ I'yr) values 
falls down to their s.p. values B(E2)w- The transition from axially symmetric to nonaxial shapes occurs due to the 
alignment of this two-quasiparticle component along the axis of the collective rotation. There the SSB effects display 
a single-particle mechanism. 
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V=1 



v = 2 





[ MeV ] 


FIG. 13: (Color online) The same as for FiglHl Two-quasiparticle energies, indicated by arrows with numbers 1, 2, 

and 3 in circles, originate from two s.p. Nilsson states (nh|[651] |[651]), (nh|[521] |[521]), and (pp|[523] ^[523]) at = 0 
and Ar = 0. 


TABLE I: The structure of the v =1 positive signature phonon in 


1 = 2 

nn|[651] 1 [660] 49% 

nn|[651] §[651] 40% 

pp|[523] ^[505] 13% 

1 = 6 

nn|[651] §[651] 72% 

nnf[651] i[660] 20% 

pp|[523] |[523] 4% 

1=10 

nn|[651] §[651] 96% 

pp|[523] |[523] 2% 

nnl[651] i[660] 1% 

I =14 

nn|[651] §[651] 85% 

nn|[521] §[521] 7% 

pp|[523] |[523] 7% 

1=18 

nn|[521] §[521] 70% 

pp|[651] §[651] 21% 

pp|[523] |[523] 5% 


In contrast with in the phonon excitation remains collective even in the backbending region 

(see Table m for / « 14?i). The B{E2; lu = 1 ^ I'yr) values decrease with the increase of the rotational frequency 
but the phonon holds the collectivity even at the transition point Klc = 0.301 MeV. Although two-quasiparticle states 
align their angular momenta along the axis x (collective rotation), the axial symmetry persists till the transition point 
(see also discussion in RefS). The mode blocks a transition to the triaxial shape. At the transition point the phonon 
becomes anomalously soft and we obtain a change of the symmetry of the HB solution from axially symmetric to 
nonaxial shape. Thus, the SSB effect occurs due to vanishing of the collective ^-vibrational excitations in the rotating 
frame. 

The second positive signature nonspurious RPA solution, ‘hiiJv =2 can be associated with the (3 - band at small 
rotational frequencies. For ^^®Dy the agreement between the experimental hug and the calculated hu ^=2 is good (see 
Figtnil. In ^®^Yb only one lowest level of /3-band is experimentally observed at hfl = 0 MeV and it is very close to 
the calculated value of the hui,= 2 - Our results may be considered as a theoretical prediction for the behaviour of the 
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TABLE II: The structure of the u = 1 positive signature phonon in 


7 = 2 

nn|[651] i[660] 48% 

pp|[523] |[523] 17% 

nni[530] |[523] 16% 

7 = 6 

nn|[651] |[651] 63% 

ppi[550] i[541] 15% 

nni[530] §[523] 13% 

7 = 10 

nn|[651] |[651] 60% 

ppi[550] |[523] 17% 

nni[530] §[523] 14% 

7 = 14 

nn|[651] §[651] 65% 

pp|[523] |[523] 18% 

nn|[530] |[523] 13% 

7= 18 

nn|[651] §[651] 71% 

pp|[523] i[541] 12% 

nni[530] §[523] 13% 


/3-band at larger rotational frequencies. 


IV. SUMMARY 

We develop a practical method based on the cranked Nilsson potential with separable residual interactions for the 
analysis of the low-lying excitations near the yrast line. In contrast to previous studies of low-lying excitations at high 
spins (cf Refs. 1^1^), we pay a special attention to the self-consistency between mean field results and description of 
low-lying excitations in the RPA. We accounted for the AN = 2 coupling in generating the Nilsson states and included 
the Galilean invariance restoring piece according to the prescription of Refl^ Moreover, we enforced the HB stability 
conditions, provided by Eq. m, that yield deformation parameters very close to the self-consistent values. Finally, 
we fixed the strength parameters of the interaction so as to ensure the separation of the spurious modes from the 
intrinsic excitations at each rotational frequency. This way we provide a reliable approach to study the spontaneous 
breaking effects of continuous symmetries of the rotating mean field. Note that even in self-consistent calculations 
with effective forces (Gogny or Skyrme type) this separation is not guaranteed due to a finite size of the configuration 
space. One needs an extended configuration space to ensure a good separation of a spurious contribution to the 
intrinsic wave function (see, for example, discussion in Ref l 2 lll . 

We analyse the rotational properties of the yrast and low-lying positive signature excitations in the transitional 
nuclei ^^®Dy and ^®^Yb undergoing the backbending. The agreement between our results and experimental data is 
remarkable. We obtain a simple expression, Eq. lEll . for the reduced transition probability along the yrast line, that 
naturally explains the increase/decrease of i3(£'2)-transitions due to the shape transition from axially symmetric to 
nonaxial shapes with different sign of the y-deformation of the rotating mean field. The magnitude of the B(E2) 
transition probability along the yrast line increases with the angular momentum and drops down at the transition 
point in both nuclei. We demonstrate a good agreement between the dynamical moment of inertia calculated at 
the mean field level and the Thouless-Valatin moment of inertia calculated in the RPA in the realistic calculations. 
This is one of the stringent test of the self-consistency between mean field and RPA calculations. Our RPA analysis 
reveals a mechanism of the backbending phenomena in rotating nuclei, which is caused by the disappearance of the 
7 -vibrations in the rotating frame and following alignment of the two-quasiparticle components of the 7 -phonon. 

We found that in axially symmetric nuclei in the ground state two types of quantum phase transitions may occur 
with rotation, which are associated with the backbending. In ^®®Dy we obtain that 7 -vibrational excitations {K = 2) 
tend to zero in the rotating frame with the increase of the rotational frequency 12 , in close agreement with experimental 
data. Although two-quasiparticle states align their angular momenta along the axis x (collective rotation), the axial 
symmetry persists, since the vibrational mode blocks a transition to the triaxial shape. Near the transition point flc 
there are two HB minima with different shapes: axially symmetric and strongly nonaxial. At the transition point the 
phonon energy tends to zero and we obtain a change of the symmetry of the HB solution from axially symmetric to 
nonaxial shape, which corresponds to the backbending in ^®®Dy. A drastic change of the mean field configuration leads 
to large fluctuations of the dynamical moment of inertia at the transition point. The observed phenomenon resembles 
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very much the structural phase transition discussed within the anharmonic Landau-type model in solid state physics 
|4fi|. We propose to consider the onset of the j-instability at the transition point associated with the backbending, 
accompanied with large fluctuations of the dynamical moment of inertia, that exceed by few times the value of the 
kinematical moment of inertia, as a manifestation of a shape-phase transition of the first order. In contrast with 
^®®Dy, at the vicinity of the transition point flc in a single neutron two-quasiparticle component dominates 

96%) in the phonon structure. There the backbending is caused by the alignment of the neutron two-quasiparticle 
configuration. The energy Eq(P,j) and the order parameter 7 (see FigEJ are smooth functions in the vicinity of 
the transition point flc. The smooth behaviour of the energy at the transition point implies a small amplitude of 
fluctuations of the dynamical moment of inertia. Extending the Landau-type theory for rotating nuclei, we proved 
that this transition can be associated with the second order quantum phase transition. Thus, if the amplitude of the 
fluctuations of the dynamical moment of inertia at the transition point is of the same order of magnitude as the value 
of the kinematical moment of inertia, a backbending may be associated with a quantum shape-phase transition of the 
second order. We expect that different mechanisms of the backbending should affect wobbling excitations that attract 
a considerable attention nowadays. This is a subject of the forthcoming paper. 
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